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INTRODUCTION 


For  a  nonequilibrium  beam-plasma  system,  the  lowest-order  Born  approxi¬ 
mation  for  the  nonlinear  force  F^,1  ^  on  the  bare  charge  of  an  unperturbed 
relativistic  test  particle  is  given  by*-3 


■f  + 

+( 1 )  2 it  (  ^  +  va*Ek  r(  +  * 

Fa  =  lim  —  ea  J  dk  k - -  k*vaj 

t  ■*■<*>  c  a)  +  i  6 


Here  a  is  a  label  designating  the  species  of  the  test  particle  of  charge  eQ 
and  velocity  va*  t  is  the  time,  k  is  the  wave  vector,  and  w  is  the  fraquency 
corresponding  to  the  Fourier  component  Ej^  of  the  total  field,  S(x)  is  the  one¬ 
dimensional  Dirac  delta  function,  5  is  a  small  imaginary  part,  and 
dk  =  d3!?  do).  Equation  (1)  differs  from  equation  (7)  of  Akopyan  and 
Tsytovich*  by  an  additional  factor  of  (2tt)_3  appearing  in  the  latter  because 
of  the  different  Fourier  transform  convention  chosen  there. 2  If  the  part  of 
the  contribution  to  the  force,  equation  (1),  which  is  associated  with  recoil 
of  the  unperturbed  bare  test  particle  due  to  bremsstrahlung,  is  averaged  over 
the  random  phase  of  the  bremsstrahlung  field,  it  reduces  to 

<Fgrad(1)>  =  <F“fad(1a)>  +  <£orad(1b)>  _  (2) 

where 

+orad(1a)  1  ?  ?  ,  d3Pfl  dw  d^  d3<  k 

<<X  ”  16*2  e«e6  1  (27T)  3 

X  lE^^^^e^Aiftk^OGipt-^v^l2  (3) 


x  (<  -  ^)*^pefp(30)6(w  -  +  i<  -  ^We) 

is  defined  and  calculated  in  this  report,  and  <F^rad(1t>)>  is  defined  by  equa¬ 
tions  (25)^  (26),  and  (47)  below  but  calculated  in  another  report.*  The 
quantity  <F^rad^  b^>  is  not  mentioned  by  Akopyan  and  Tsytovich,  1  but  it  is 
nonvanishing  and  makes  an  important  contribution  to  nonlinear  bremsstrahlung. 


*A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plazmy,  _2_  (1975),  673  [Sov.  J.  Plasma  Phys.  (1975),  371], 

2ff.  E.  Brandt,  Nonlinear  Force  on  an  Unpolarized  Relativistic  Test  Parti¬ 
cle  to  Second-Order  in  the  Total  Field  in  a  Nonequilibrium  Beam-Plasma  System, 
Harry  Diamond  Laboratories ,  HDL-TR-1995  f August  1983). 

*H.  E.  Brandt,  Nonlinear  Dynamic  Polarization  Force  on  a  Relativistic  Test 
Particle  in  a  Nonequilibrium  Beam-Plasma  System,  Harry  Diamond  Laboratories , 
HDL-TR- 1994  (September  1983). 

*H.  E.  Brandt,  Collective  Bremsstrahlung  Recoil  Force  on  the  Bare  Charge 
of  an  Unperturbed  Test  Particle  in  a  Nonequilibrium  Beam-Plasma  System,  Harry 
Diamond  Laboratories ,  1983,  (to  be  published) . 
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The  quantity  <F2ra<^  a^>  given  by  equation  (3)  corresponds  to  equation  (27)  of 

Akopyan  and  Tsytovich. 1  Their  equation  (27)  is  expressed  in  Gaussian  units, 
whereas  equation  (3)  above  is  in  MKS  units.  This  results  in  an  additional 
overall  factor  of  (4it)2  in  Akopyan  and  Tsytovich.1  Also  there  is  an  addition¬ 
al  factor  of  (2it)“6  in  their  work  because  of  different  Fourier  transform 
conventions.  However,  their  equation  (27)  also  differs  from  equation  (3)  here 
in  that  they  omit  an  overall  factor  of  2.  Also,  the  arguments  k  and  -tc  of 
(k,-x)  in  equation  (3)  appear  incorrectly  interchanged  in  equation  (27) 
of  Akopyan  and  Tsytovich. 1  As  shown  below,  these  disparities  are  evidently 
due  to  minor  errors  in  their  work. 


In  equation  (3),  0  designates  the  species  or  tne  otner  particle  witn 

velocity  vg,  momentum  pg,  and  background  distribution  fpl°^  involved  in  the 
bremsstrahlung  process.  Also,  ic  is  the  momentum  transfer  from  the  test  parti¬ 
cle  and  E®^)  is  the  amplitude  of  the  lowest  order  stochastic  bremsstrahlung 
field  with  polarization  e^  and  mode  o.  The  matrix  Af^(k,-<)  in  equation  (3) 
is  defined  by1”3 


(ki,k)  = 


5i]  + 


v°iV^rki  -  jv 

/  -►  "  \n 

-  k  •Vq  -  iS  (w  -  k‘va  -  i6J * 


The  quantity  Gmn(k)  in  equation  (3)  is  the  linear  photon  Green's  function, 
whose  inverse  is  given  by4 


(kj.kj  -  k2«ij  ) 


(5) 


1A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plazmy,  1_  (1975),  673  [Sov.  J.  Plasma  Phys.  1_  (1975),  371]. 

2//.  E.  Brandt,  Nonlinear  Force  on  an  Unpolarized  Relativistic  Test  Parti¬ 
cle  to  Second-Order  in  the  Total  Field  in  a  Nonequilibrium  Beam-Plasma  System, 
Harry  Diamond  Laboratories ,  HDL-TR-1995  (August  1983). 

3H.  E.  Brandt,  Nonlinear  Dynamic  Polarization  Force  on  a  Relativistic  Test 
Particle  in  a  Nonequilibrium  Beam-Plasma  System,  Harry  Diamond  Laboratories , 
HDL-TR- 1994  (September  1983). 

E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi¬ 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories ,  HDL-TR- 
1996  (September  1983). 
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For  a  spatially  isotropic  system,  the  dielectric  permittivity  tensor  e-  •  is 
given  by4 *”7  13 


+  /  k  j_k -»  1  ->•  kik-;  ^ 

ei  j  (k,  oj)  =  ^63.  j  ~  '~^2~ t^k,ti))  +  Ejj(k,(o)  ,  (6) 

where  et  and  are  the  transverse  and  longitudinal  permittivity,  respec¬ 
tively.  In  this  case  the  linear  photon  Green's  function  is  given  by4~6 7 


«ij  - ^ - j(«ij  -  , 

k2e£  Et ( w  +  i6)2  -  ^  \  ^  / 

and  the  transverse  and  longitudinal  permittivities  are  given  by 
1  1  v  es  r  d3Ps  ,  k2V2  -  ajt*vs 

£t  -  e° +  2  77777)2  3  72  ^  TTo3  i  Vsm  +  777X777 km_ 


_J _  y  fl  r  llPs _ 1_ _ 

“  +  15  s  (2tt)3  u)  -  )t.vs  +  i6 


R(0) 

Ps 


(9) 


respectively.  Here  eQ  and  uQ  are  the  permittivity  and  permeability,  respec¬ 
tively,  of  the  vacuum.  Each  species  s  of  particles  present  in  the  system 
makes  a  contribution  to  equations  (8)  and  (9).  If  the  distribution  function 
depends  only  on  energy,  then  using  the  expression  eg  for  the  relativistic 
energy  of  a  particle  of  species  s,  namely. 


es  =  (Ps°2  +  role4)1/2  ,  (10) 

allows  equations  (8)  and  (9)  to  be  easily  rewritten  in  terms  of 

3^0)/3es.4.5 


4fl.  E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi¬ 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories,  HDL-TR- 
1996  (September  1983). 

-*A.  J.  Akhiezer ,  I.  A.  Akhiezer ,  R.  V.  Polovin,  A.  G.  Sitenko,  and  K.  N. 
Stepanov,  Plasma  Electrodynamics,  Vol .  1,  Linear  Theory,  Pergamon  Press 
(1975),  206. 

*>V .  N.  Tsytovich,  Nonlinear  Effects  in  Plasma,  Plenum  Press,  New  York 
(1970),  314. 

7V.  N.  Tsytovich,  Theory  of  Turbulent  Plasma  (Consultants  Bureau),  Plenum 
Publishing  Corp .,  New  York  (1977),  63-65. 
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In  the  present  work,  equations  (2)  and  (3)  are  derived  from  first  princi¬ 
ples  for  the  purpose  of  identifying  some  of  the  important  physical  assumptions 
and  mathematical  techniques  involved  in  Tsytovich's  theory  of  nonlinear  brems- 
strahlung  and  radiative  instability  in  relativistic  nonequilibrium  beam-plasma 
systems  (see  also  Selected  Bibliography).1  This  will  facilitate  understanding 
and  illustrate  the  procedures  for  reducing  other  components  of  the  total 
nonlinear  recoil  force  on  a  test  particle  due  to  bremsstrahlung  emission.  1 

The  present  work,  together  with  other  related  work  by  the  author  (see  also 
Selected  Bibliography ), 2  4  is  important  for  ongoing  work  in  calculating  col¬ 
lective  radiative  processes  and  conditions  for  the  occurrence  of  radiative 
instability  in  relativistic  nonequilibrium  beam-plasma  systems. 

2.  REDUCTION  OF  A  COMPONENT  OF  THE  BREMSSTRAHLUNG  RECOIL  FORCE 

In  this  section,  equation  (3)  is  to  be  derived  by  reducing  that  part  of 
equation  (1)  which  corresponds  to  bremsstrahlung  recoil.  Also  the 
quantity  >  appearing  in  equation  (2)  will  be  defined.  It  has  been 
shown  by  using  the  balance  equations  that  the  general  form  for  the  collective 
bremsstrahlung  recoil  force  on  a  test  particle  participating  in  induced  brems¬ 
strahlung  in  a  nonequilibrium  beam-plasma  system  is  given  by1'8 

■  -4 * * *’3  /  -  *)-(Vp6)  ("’ 


a)2  3w 


X  |e£(°>|2  v£arP0(*,J?><s(cu  -  if.*a  t  (if  - 


Ri 
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E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi¬ 
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1996  (September  1983). 
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tories  HDL-TR-2026  (January  1984). 


a;av.vv.v>;.' \-\v; \  -  v  *  \  - 


Here  the  dielectric  constant  e°  (k,co)  for  the  bremsstrahlung  mode  a  is  defined 
in  terms  of  the  dielectric  permittivity  tensor  e,j(k,u))  and  the  unit  electric 
polarization  vector  <££  by4'6'8'* 


e  (k,  a)) 


a*  a  c2  f*  +a\r+  +o*\ 

ekieijekj  +  eo  lk,ekKk,ek  J  • 


(12) 


Also  in  equation  (11),  Ejj^O)  is  the  zeroth-order  amplitude  of  the  stochastic 
bremsstrahlung  field,4  and  Vp  fpg(£»£)  is  the  bremsstrahlung  transition  proba¬ 
bility  with  the  delta  function  factored  out,  expressing  conservation  of 
energy.  It  is  clear  that  vpa,po^K*^)  must  be  at  least  first  order  in  the 
regular  part  of  the  field.  In  reducing  equation  (1),  only  those  contributions 
which  are  of  the  same  form  as  equation  (11)  are  to  be  included,  since  only 
they  correspond  to  bremsstrahlung. 

The  total  field  E^  appearing  in  equation  (1)  and  associated  with  the 
bremsstrahlung  process  involving  the  test  particle  is  given  by1'4 


(13) 


Here  e£(0) 


is  the  zeroth-order  bremsstrahlung  field 


e£<0)  =  e£(0)££ 


(14) 


where  o  designates  the  mode  and  is  the  polarization  vector.  If  the  brems¬ 
strahlung  field  is  to  be  real,  then 


e£(0)  =  E?i0)* 


(15) 
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The  fields  Ej^1  ^  and  in  equation  (13)  are  the  self  fields  arising  directly 

from  the  relativistic  test  particle's  cwn  motion,  namely,4 


(1) 

Ekn  = 


(  2tt)  3(  a)  +  i6) 


vconGnm<kl  «(“  "  k*va) 


Ekn  '  3  Gnm^k  1 

(  2tt)  J(  (d  +  l  6) 


dk  i  *  (a)*,-  i  +  +  -. 

x  J  ^~T~n  Eki  j  Amj  lk1'kJ6l“  +  a.,  -  k.va  -  krvaj  . 

The  fields  and  E^  in  equation  (13)  are  the  increasing  order  fields 

produced  by  the  dynamic  polarization  current  induced  by  the  test  particle  and 
are  given  by4 

( 1 )  i  dk!  dko  , 

Edpkn  “  Gnm<kl  1  es  /  ,  .  *w  '  -i  +  kl  +  k2 ) 

2  (  u  +  i  6 )  s  (uh-i6J[id2-i6J 


x  [smj  Jl(k' "k1  »  ~k2  )  +  smH]  (k'  -k2»  _k1  )  ]Ek-j  jEk2£ 


Edpkn  “  Gnm(kl  1  e si 

0)  +  16 


dk-j  dk2  dk3 

(o>1  -  i  6  ]  ( u>2  -  i 6 ) ( 003  -  iS) 


x  6(k  +  k!  +  k2  +  k3)  jj,(k,  -k  i  f -k2 ' -k3 )E&-|  iEic2 jEk3 i.  * 

respectively.  The  symmetrization  in  equation  (19)  appearing  in  earlier 
work1'4  has  been  removed.  The  delta  functions  appearing  in  equations  (18)  and 
(19)  are  four-dimensional,  namely, 

S(k)  =  63(ft)6(u>)  .  (20) 


The  second-order  nonlinear  conductivity  tensor  s|  j  ]j(k»ki /k2  )  appearing  in 
equation  (18)  is  given  by 

1A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plazmy,  1_  (197 5),  673  [ Sov .  J.  Plasma  Phys.  1_  (1975),  371]. 

4H.  E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi¬ 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories,  HDL-TR- 
1996  (September  1983). 


■VVV1* 
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^  £> 


’Si 


sijjl(k,k1,k2j  =  e|  /  -  — -  - - —  Oi  -  k^Vg) 

(2ir)  u>-k*vs  +  io[_ 


+  v-ik 


3*1m  3p« 


m]  [isl 


vsl  3 

9ps 


A  • 

1  R(0) 

Jfps  • 


The  third-order  nonlinear  conductivity  tensor  £i}lm(k'k1 ,k2»k3)  is  given  by 
*•!  jlm(k,k1  'k2»k3) 


=  -iel  / 


( 2tt  ) 3  a)  -  k*Vq  +  i6 


x  [6jn(<*>1  -  ki*vs)  +  k1nvsj]  -r| - 7?1— -»  T  y - ~ 

dPsn  t*J  —  oj i  -  [k  -  kij»vs  +  16 

x  [6fu(“2  “  £2*vs)  +  k2uvsl]  -A - V"  T - — 

dPsu  0) 3  -  k3*vs  +  16 


x  [6mqC“3  "  k3*vs)  +  k3qVsm] 


i  3  ,R( 0) 

*lVs®J  3psq  fPs 


Letting  BF^,1  ^  denote  the  operation  of  including  only  that  part  of  F^1  ^ 
which  is  of  the  bremsstrahlung  form  given  by  equation  (11),  we  obtain 

<Jarad(1)>  =  0<|H  )>  .  (23) 

Here  <  >  denotes  the  ensemble  average  over  the  statistical  phase  distribution 
of  the  bremsstrahlung,  which  to  the  needed  order  is  assumed  to  be  random.*4 
Substituting  equations  (1)  and  (13)  in  equation  (23)  gives 


<fcrad(1)>  =  J  ^arad(1)>n  f 
n=1 


where 


iorad  ( 1  )  at  r  an  *  ■*  ♦n  . ,  ♦  ♦  •> 

<Fa  >n  =  B  T~  ea  /  - 7  k  va»<Ek>6(u  -  k«va)  (25) 

t+*  *-  uj  +  io 

**//.  E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi¬ 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories ,  HDL-TR- 
1996  (September  1983). 


dk  ♦  -*■  +n 


and 


{Sg,  n  =  1,  2,  ...5}  =  {e£(0),  W>,  E^2),  E^,  . 


To  the  needed  order,  the  stochastic  properties  of  the  bremsstrahlung  field  are 
approximated  by 1 ' 4 


<e£)01>  =  0 


-  •&«&*  k(O,i2«0‘  *  *i) 


Substituting  equations  (26)  and  (27)  in  equation  (25)  for  n  =  1  produces 


<F^rad( 1 )>i  .  0  . 


If  we  substitute  equations  (26)  and  (16)  in  equation  (25)  for  n  =  2,  it  is 
evident  that  the  resulting  expression  is  not  of  the  required  form  ^,iven  by 
equation  (11).  For  example  it  does  not  contain  |e£  °*|2»  therefore  the  B 
operation  gives  zero,  thus: 


£orad(  1 ) 


>2  =  0 


Next  substituting  equations  (26)  and  (17)  in  equation  (25)  for  n 
obtain 


=  3,  we 


►orad ( 1 ) 


>3  =  B  lim 
t-*» 


2tt  r  dk  >  \ 

ea  j  kvan6(a)  -  k*vaJ  — y- — — ~ 

u)+i6  [_(2tt)3(uj  +  io)_ 


- 77  <Eki  j >^mj  [k, ,kJ6(a)  +  a),  -  k«va  -  k^vj 

a>i  -  16  J 


*A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz,  Plazmy,  1_  (1975),  673  l Sov .  J.  Plasma  Phys .  1_  (1975),  371). 

E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi¬ 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories,  HDL-TR- 
1996  (September  1983). 
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Furthermore,  substituting  equation  (13)  in  equation  (31)  and  using  equation 
(26)  then 


iorad(l)  _  y  iorad(l) 

<Fa  >3  -  L  <Fa  >3n  ' 

n=1 


where 


<F 


-*arad ( 1 ) 


„  ,  .  2n  t  dfc  *  ,(  *  i 

>3n  =  B  lim  —  ea  /  - —  kvotn6[a)  -  k*vaJ 

t+»  1  to  +  1 5 


dk 


a 


_e  a^nm  ^  ^  ^ 


(  2tt)  3(  to  +  i  6) 


/ 


dk! 


(32) 


(33) 


n*  (a)*/  \  ✓  \ 

.  <Ek1j>Amj  [k1fkJ6(io  +  to!  -  k*va  -  k!  *vaJ  . 
0)^  -  i  o  1 

Substituting  equations  (26)  and  (27)  in  equation  (33)  for  n  =  1  produces 


<^orad(1)>31  _  o  , 


(34) 


Next  substituting  equations  (26)  and  (16)  in  equation  (33)  for  n  =  2,  we 
obtain  an  expression  which  is  clearly  not  of  the  bremsstrahlung  form  given  by 
equation  (11),  and  therefore  the  B  operation  yields  zero,  thus: 

<Fgrad(1)>32  =  o  .  (35) 


Next  substituting  equations  (26)  and  (17)  in  equation  (33)  for  n  =  3,  we 
obtain  a  higher  order  Born  term  which  is  also  not  of  the  bremsstrahlung  form 
given  by  equation  (11),  and  therefore 

<#grad{1)>33  =  0  .  (36) 

For  consistency  with  the  Born  approximation,  only  those  terms  which  are  first 
or  second  order  in  the  regular  part  of  the  field  are  needed. 1 

So  that  in  the  limit  of  t  +  *  the  integrals  in  equation  (25)  or  equation 
(33)  do  not  vanish,  it  is  important  to  note  here  that  the  following 
combination  of  delta  functions  is  needed  in  the  integrand: 


$2(u>  -  k»va)  =  lim  —  6(u>  -  k»va) 


t-M» 


2tt 


(37) 


In  this  way  the  t-1  factor  is  canceled.  Equation  (37)  is  obtained  as  follows. 
The  well-known  integral  representation  of  the  delta  function  is 


6{u>) 


-iujt 

2ir 


(38) 


*A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma ,  Fiz.  Plazmy,  1  (1975),  673  [ Sov .  J.  Plasma  Phys.  1  (1975),  371]. 


Therefore 


-  k*vj  =  lim  I 

J-t 


dt '  -i  (o)-k  *v 


Then  using  equation  (39)  one  has 


*2f  t  *  i  xf  £  *  )  ft/2  dt'  -i(“-k’va)t' 

6z[a)  -  k*vaJ  =  6[a)  -  k*vaJ  lim  /  — —  e 

t+<°  J  -t/2  211 

ft/2  dt'  ,  ♦  -  x  -i(ahk.va)t' 

'  J-t/ 2  2ir  "  k*va)e 

.  ft/2  dt1  /  +•  +  s 

'  t-  it/2  (“  ‘  Va) 


6(a)  -  k»va) 


=  lim  6(u  -  k*va) 


Thus  in  evaluating  equation  (31)  one  has  an  integrand  containing  the  factor 
F  =  —  6[o)  -  k«vaJ<EjCl  j>6[oi  +  ag  -  k*va  -  ki  *vaJ 
or 

F  =  ~  6(b)  -  ^•va)<E)J1  j>6(a)1  -  ki*va)  ;  (' 

or  substituting  equation  (13)  in  equation  (42)  produces 


F  =1  T  6(a)  -  k*va)<E{J*j>6(o)1  -  ki«va) 
n=1  ' 


Because  of  equation  (37),  another  factor  of  6(u)1  -  £*va)  is  needed  to  give  a 
nonvanishing  limit.  For  n  =  2  in  equation  (43),  one  picks  up  a  factor  of 
6(a)-,  -  l^i  *va)  by  equation  (16);  however,  the  remaining  integrand  does  not 

contain  the  factor  |e^°^|2  necessary  to  be  of  the  bremsstrahlung  form  given 
by  equation  (11).  An  iteration  of  equation  (17)  in  equation  (43)  can  give  the 
quadratic  delta  function,  but  the  remaining  integrand  is  higher  order.  By 
analogous  reasoning 

<?grad(1)>34  =  0  >  (44) 

<^radm>35  =  0  .  (45) 

Therefore,  substituting  equations  (34)  to  (36),  (44),  and  (45)  in  equation 

(32)  produces 


f  *  ,  ^*4  f  #  f  ,  r  ^  •  TV*  .  *«  #  *  #  .  *  •  9  •  *  .  •  «  *  •  *  «  *  «/  •  ’  •  *  *  .*  4^  <  ’  ^  -  1  *  *  4^4*  •*»*•*  *  *  »  ’  ‘  , 


.  *  .  "  .  '  4  * 

*r. 


<?Srad(1)>3  =  0 


(46) 


The  quantity  <?^rad^1^>4  is  nonvanishing  and  is  calculated  in  a  separate 
report.*  Here  the  result  of  that  calculation  will  be  designated  <?^rad^'b^>. 
Thus 


<Fgrad(1)>4  =  <F°rad( 1b) >  . 


(47) 


In  the  present  report  we  emphasize  the  calculation  of  <Forab^1a^>  in  equation 

(2). 


Substituting  equations  (26)  and  (19)  in  equation  (25)  one  has 


<F 


♦orad ( 1 ) 


„  2ir  r  ax  \ 

>5  =  B  lim  —  ea  J  -  k6[aj  -  k*vaJv, 

t  -**>  a)  +  i  6 


dk  , 


a 


an 


-i  ~  r  f  dk1  ^2  dk3 

*  .  .  Gnm^^  1  es  /  (  w  rw  ,  i 

ai+i6  s  (ii)1-i6j[a)2-i6J((»)3-i6j 


(48) 


x  6(k  +  k^  +  k2  +  k3  )^mi]  i(k' -k1 » -k2' “k3  )<Ek-|  iEk2  jEk3£> 


One  has,  using  equation  (13)  to  the  required  order  in  equation  (48), 


<Ek1iEk2jEk3£>  =  <(Ek1(?)  +  4]i  ♦  4 +  4p)li +  ESi) 


*  (<?’  +  4'2]  *  *  4^2j  * 

-  izk'J.’  *  4\i  *  4|l  *  4y3>.  *  4p4i)>  • 


(49) 


If  equation  (49)  is  substituted  in  equation  (48),  the  only  possible  nonvan¬ 
ishing  contributions  to  the  required  order  are  given  by 


*H .  E.  Brandt,  Collective  Bremsstrahlung  Recoil  Force  on  the  Bare  Charge 
of  an  Unperturbed  Test  Particle  in  a  Nonegui librium  Beam-Plasma  System,  Harry 
Diamond  Laboratories ,  1983  (to  be  published  as  HDL  technical  report). 
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:•  ,% 


.*  «• 


£orad( 1  ) 
<-r  a 


>52 


B  Jl:  r 1  ~Te  fst“ '  k*M'«)Vop 

*(-fr 1-/7 - d^7^rr — ^ 

\0)  +  16/  S  (u>i  -  1 6  J  ( u)2  -  i  6 )  ( w3  -  i  6 ) 

x  6(k  +  k,  +  k2  +  k3  )4i  j  $(k,  -k-j , -k2,  -k3 )  1  E^j*  <e£J  S’’ *£*3°’ *> 


and 


<?Srad"  >>53 


-  B  lim  —  ea  /  - —  k6(aj  -  k'VaJvcp 

t+»  z  a)  +  i6  ^ 

\(jj  +  i6/ 


|Gpm<k>  I  es  1 


dk-)  dk2  dk3 


(a>]  -  i  6  )  ( w2  -  i  6 )  ( io3  -  i6) 

x  6(k  +  kl  +  k2  +  k3)4f]£(k,-k1<-k2,-k3)  j  4]V<*Z{2V**ZT>  • 

Substituting  equations  (28)  and  (16)  in  equation  (52)  produces 

<?Srad(1,>51 


=  B 


.  .  2ir  ,  dk  , 

lim  —  ea  /  - —  k6(u)  -  k *va) 

t+«  c  ai  +  16 


x  v 


C  *\«)Gp" 


<k>  I  es  / 
s 


dk^  dk2  dk3 


(ail  -  i  6  )  (  w2  -  i  5 )  ( iu3  -  i  6 ) 
x  6(k  +  kl  +  k2  +  k3)4fj£(k,-k1r-k2,-k3) 

X  *  [(2^1^-  i«)]VanG^l(lt3)6(U>3  '  ^3*^a)ek*iek1  j 

x  6(k-|  +  k2)  . 


a(0) 


Integrating  equation  (55)  first  over  k2,  using  the  property  of  the 
function  6(k^  +  k2)  and  effectively  setting  k2  =  -k^  (that  is,  £2  =  -Jfci 
-u1 ),  we  obtain 
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(53; 


(54) 


(55) 


delta 


1.1  I. 


;orad ( 1 ) 


>5T  =  B  lim 


-e|  dk  dk-|  dk3  k6(k  +  k 3  ) 6 ( to  -  k*va) 


8ir2t  (o)+i6)2(u)i  -  +  i6)(io3-i6)2 


x  vapGpm  <k>  I  eS^^(k.-k1,k1,-k3] 


x  '  ^3*va)e^*ie,°1  j  |e)°i(0)  j  2 


Next  integrating  equation  (56)  over  k^  produces 


+orad( 1 ) 
a 


-e&  ,  dk  dk  i  +  *  + 

>51  =  B  —  o-  /  r~~~ir4~2  k(S  -  k#va) 

t-*-oo  87r^t  (OJ  +  1<5)HC0^ 

X  vapGpm  (k)  £  es  ^mi]  i(k»  “k  1 1  k  1  ,  k  jv^G  -k  ) 


e,a*  e,°  ■  eJ{0)  2 

fcki 


One  also  has 


Gnm  ^ -k  ^  Gnm ( k )  . 


Equation  (58)  holds  since4 


Ekn  ~  .  .  Gnm(k)j)cm 

U)  +  10 


and  by  the  reality  of  the  current  and  field 


E-kn  =  Ekn 


j-kn  =  3kn 


E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories,  HDL-TR 
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r.  -VV.  • 
r.V>  A  .>  V. 


•/  *  *  •/  «/  «/ v  *  *  .  - .  ^  ■  .  »  .  <  .  » ,  • .  * .  «  , «  ,  A  ^  /,  .  .  ,•  .»  ,• 


Talking  the  complex  conjugate  of  equation  (59),  replacing  k  by  -k,  and  substi¬ 
tuting  equations  (60)  and  (61)  produces 


®kn  ~  .  Gnm^'k^Okm  *  (62) 

to  +  16 

Comparing  equations  (59)  and  (62)  shows  that  equation  (58)  follows.  Next 
substituting  equations  (37)  and  (58)  in  equation  (57),  and  noting  that  the  t's 
cancel,  we  find  that  the  limit  is  trivial  and 


>orad(  1 )  /~ea  \  f  <*ki  kfi(o)  -  k*va) 

<F“  >5'  ’  %Ti)  1  - 

x  vapGpm()c)  I  es  ^mij  i  0* •  1  1  *  ^  )vc<n^£ji  )  (63) 

s 

Changing  variables  {k,k1  }  to  jk^k)  and  |i,j,£}  to  { £, i ,  j  }  in  equation  (63) 
produces 


♦arad  ( 1 )  (~e  a  \  r 

<p“  >5'  *  B\TT,V  I 


!  dk  dki  |  „o( 0)  |  9 

-si  o ,  .  k1  lEk  2 


+  ifi)1* 


x  vapGpm(k1  JvcGjnfk!  JejfJegj. 
x  Jti j  (k1  »_k,^,ki  )<5(a)i  -  k-|  'Vqi) 

If  we  substitute  equations  (28)  and  (16)  in  equation  (53),  we  obtain 

<?»~d'">52 

=  B  lim  ea  /  — — —  k6(oo  -  k*va] 
t  a)  +  i6 

/  -i  V  .  v  (  dk1  dk2  dk3 

*  Va%  ♦  16/  pm  k  6s  J  («,  -  i  6 ) ( uj2  -  i«K»3  -  IT) 

x  6(k  +  k,  +  k2  +  ^3)4ijA(k»-k1,-k2,-k3)  \  jj~ ~ ~ j 
x  vajGjn(k2)6(<^2  “  ^2  •va)e]?*ie£1  «|  Efc{  0 )  |  ^(k,  +  k3)  . 


If  we  integrate  over  k^  and  k2»  equation  (65)  becomes 


♦or ad ( 1  )  „  , . 

<F a  >52  =  B  lim 

t  ■+•» 


(  ~e0L  \  r  ^  <^1  k 

im  l  — —  I  J  - — 

-►oo  \8Tt2t/  (to  +  i6)4u 


x  62(uj  -  k  •vot)vQpGpm(k )  £  es  E^ii  £, (k ,  -k^  ,k,  ki  } 


x  vonGjn(-k)ekiieki £  Ek^°^  2 


Next  substituting  equations  (37)  and  (58)  in  equation  (66)  produces 


♦orad  ( 1 ). 


/~ea  \  r  ^  dk1  k  ,/  >  -►  . 

■BW)  1  FTTTjx  i(“-k‘v“) 

x  vapGpm(k)  I  es  ^mij  £  (k '  -k  1 •  k'  k1  )vanGjn  ^k  J 

x<iek°1£|Kka1(0)|2  • 


Changing  variables  {k.j,k}  to  {k , k ^  }  and  { i , £ }  to  { £ , i }  in  equation  (67),  we 
obtain 


+orad(1 ) 


x  vopGpm(k1  )vonGjn(k1  )ek£eki 
x  I  eszm£ji(k1'-k'k1'k)5(a>1  "  *WcJ  • 


Next,  substituting  equations  (28)  and  (16)  in  equation  (54),  we  get 


♦arad ( 1 ) 

<Fa  >53 


dk  ,♦  * 


B  lim  ea  /  k 6(a)  -  k-vj 

t+o°  111  +  10 

(  -i  \  v  f  dk-j  dk2  dk3 

*  v,p\“r«jGp”(k)  i  *s  J  'i^~n)(o2  -  n)^ 

*  S(k  *  k,  +  k2  *  k 3 ) j  , (k ,  -k i ,  ~k 2 1 -k 3 ) 

-  5  [(»,)“(!!,  -  )5(“’  - ki 

»=k°;jek°2iiikki0li2«(i'2  +  k3)  . 


mmm . 


>  »  »  . 


(70) 


If  we  integrate  over  k^  and  k1 ,  equation  (69)  becomes 


♦or ad ( 1 ) 

<F„  >53  .  B 


u„  (--4)  I 

t+®  \8ir2t/  (a)  +  i6)4oj 


x  52((0  -  it-VaJv^GpnOO  |  «s EiSii j »(*•*■  -*2'k2 ) 
-v1„G*n(-1Oe1«;j.g2l|Ek«lt»|2  . 

Next  substituting  equations  (37)  and  (58)  in  equation  (70)  produces 


>orad  ( 1 )  „/~ea  \  r  d 

<F«  >53  ■  BW) 1  J7:i 


dk2  k  v 

77771  «(“-  1“v“) 


x  vapGpm  (k)  I  es ^mi]  i(k' k» -k 2> ^2 ) 


x  vanGin(k)ek2jek2£  Ek2°^  2 


Changing  variables  {k,k2}  to  {k  » k }  and  {i ,  j  ,  A }  to  { j  ,  £,  i  }  in  equation  (71) 
produces 


.*orad(1).  „  (~e  a  \  ^  dk,  |„o(0>| 

<F“  >53  -  B(^sJ  1  Zh^TUW  1  k  1 


”  Vl6,3/  '  „2(«,,  ♦  16)-  1  *  1 

x  VapGpm  (k,  jVanGjn  (k,  je^ej^  ( 72 ) 

x  I  «sJmjli(k1»k1'-k'k)6(®1  • 

s 

However,  equation  (72)  is  not  of  the  bremsstrahlung  form  given  by  equation 
(11);  in  particular,  it  does  not  contain  the  delta  function  factor 
6(w  -  +  ($  “  in  the  integrand.  Such  a  factor  could  only  appear 

through  the  third-order  nonlinear  conductivity  tensor  appearing  in  equation 
(72).  According  to  equation  (22),  the  factor[u>  —  uj^  -  (1?  -  l?i)»\Fs  +  i6]_1 
appears  in  Eli|m(k,ki  ,k2,k3),  and  it  can  be  written 


u  -  a>i  -  (k  -  k)  )»vs  +  i<5  u-(Di  -(k-k))' 


in6(oo  -  o>i  -  [J?  -  j?i  )*vs)  , 


,  (.  •  . •••.  ■•••.*  *  *  S  '.  -  .  •  .-/•  *>••  '•  ' 

«  *  .  •  •  *  •  m  •  *  -  _•  ,  -  .  .  .  <  _  «  *!.«.»  , 


^••**'k**"«  '  I  •  m  9  »  *  ,  ’  »  Vh  *  » 

*  N’.V'.VV/  /.V.v'.  -  . 

•  ,  «  ,  v,  ,  »  „  - 


and  the  imaginary  part  contains  a  delta  function;  however,  in  equation  (72)  it 
is  ^  (ki  ,k-|  ,-k,k )  which  appears,  so  that  the  delta  function  becomes 
— i it 6 { 0 ) ,  which  is  not  of  the  required  form.  Thus  the  B  operation  in  equation 
(72)  yields  zero,  and  therefore 


.jjorad  ( 1  ) 
<Fa 


>53  =  0 


(74) 


Next  substituting  equations  (64),  (68),  and  (74)  in  equation  (51)  and  the 

result  together  with  equations  (29),  (30),  (46),  and  (47)  in  equation  (24), 

and  comparing  with  equation  (2),  then 


orad  (la) 
<Fa  > 


\  y  e  ,  ^k1  ^1 

^1 6  7r  y  s  S  oj2  ( a)i  +  i  6 ) 4 


2 


x  vapGpm0l  )vanGjn(k1  )ekieki6(w1  "  *1  *va) 


(75) 


x  [^m£ij  0*1  )  +  ^m)l]i  (^1 '“^'^1 )  1  * 


It  remains  to  be  shown  in  this  report  that  equation  (75)  reduces  to  equation 
(3). 


Evidently  equation  (75)  should  agree  with  equation  (24)  of  Akopyan  and 
Tsytovich. 1  Clearly  there  is  a  typographical  error  in  the  latter  since  a 
factor  of  ea  has  been  omitted.  Also  a  factor  of  (-2)  and  complex  conjugation 
signs  on  the  conductivity  tensor  have  apparently  been  omitted  as  shown  below. 
Their  equation  (24) 1  contains  an  additional  overall  factor  of  (4x)2(2it)_6  due 
to  the  use  there  of  Gaussian  instead  of  MKS  units  and  differing  Fourier 
transform  convention. 


As  was  already  noted  elsewhere,4  extra  factors  of  (4n),  (2tt)“6,  and  (2n)~3 
evidently  enter  through  in  equation  (19),  due  respectively  to  the  use  of 

Gaussian  units,  differing  Fourier  transform  convention,  and  differing  normal¬ 
ization  of  the  background  distribution  fp^0^  in  equation  (21).  However  the 
normalization  in  equation  (24)  of  Akopyan  and  Tsytovich1  is  apparently  the 
same  as  that  here.  Otherwise  an  additional  factor  of  (2ir)~3  was  omitted  there 
in  equation  (24).  A  factor  of  (2n)3  arises  if  the  normalization  is 
restored.  An  additional  factor  of  (4tt)(2tt)3  enters  through  in  equation 

(16)  due  to  the  use  of  Gaussian  units  and  differing  Fourier  transform 
convention.  As  already  noted,  an  additional  factor  of  (2u)-3  enters  through 
equation  (1)  due  to  the  differing  Fourier  transform  convention.  These  factors 
combine  to 


[  (4n)(2n)"6(2if)“3][(2n)3][4n(2ir)3][  (2ir)'3]  =  (4ir  )2  (2ir  )"6  . 
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Because  uie  same  definition  of  E£x  in  equation  (28)  here  is  used  in  equation 
(13)  there,  the  differing  Fourier  transform  convention  is  compensated  for. 
Thus  equation  (75)  is  in  essential  agreement  with  equation  (24)  of  Akopyan  and 
Tsytovich1  up  to  the  omission  of  a  factor  of  (-2ea)  and  a  complex  conjugation 
sign  on  the  conductivity  tensor.  This  follows  since  the  6  in  +  i<$  may 
clearly  be  ignored,  and 


Also  if  the  background  is  assumed  to  be  nondissipati  ve  to  the  needed  order 
then  it  follows  that 

((0!  +  i6)-1Gnj(k1)  =  [(«o,  +  i6)-lGjn(k1)]*  ,  (77) 

which  is  true  under  certain  conditions  to  be  discussed  below.  Also  for  a 
spatially  isotropic  system  the  Green's  function  is  symmetric  in  its  indices 
according  to  equation  (7).  It  should,  however,  be  noted  by  comparing  equation 
(14)  of  Akopyan  and  Tsytovich1  with  equation  (21)  of  Brandt1*  that  the  indices 
of  G^j  are  interchanged  in  Akopyan  and  Tsytovich1  in  the  definition  of  the 
Green's  function.  It  also  follows  that  if  one  ignores  the  zero  frequency 
resonance  denominator,  then  equation  (77)  may  be  used  provided  the 
contribution  of  the  imaginary  part  of  the  photonic  Green's  function  is 
negligible  to  the  needed  order.  Evidently  the  complex  conjugation  on  the 
conductivity  in  equation  (24)  of  Akopyan  and  Tsytovich1  has  been  omitted.  It 
ls  required  if  their  equation  (25)  is  to  be  consistent  with  the  fact  that  the 
first  and  second  complex  denominators  in  their  equation  (21)  are  (w  -  (£*vs) 
-  i 6 )  and  (w  +  oi1  -  (it  +  l?i)*vs  -  i 6 1 ,  respectively.4  In  equation  (76)  the 
tensor  Ti jlm(k'k1 »k2'k3 )  is  given  by1'4 
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It  was  concluded  elsewhere4  that  for  consistency,  negative  imaginary  parts 
must  be  implicit  in  the  first  two  denominators  of  equation  (21)  of  Akopyan  and 
Tsytovich. 1  By  equation  (22)  here  it  is  evident  that 


4flm(k1 .*,-k,-ki )  •  (79) 

Also  from  equations  (22)  and  (78)  it  follows  that 

4]lm(_k1 ,k'  -k'  -k1  )  =  “ir^i]lm(k1 ,k' _k' ~k1  )  *  (80) 

If  we  substitute  equation  (80)  in  equation  (79),  then  equation  (76)  follows. 

Equation  (77)  holds  approximately  if  there  is  negligible  linear  absorption 
of  energy  by  the  background  beam-plasma.  To  see  this  one  notes  that  the 
energy  ej[bs  absorbed  linearly  by  the  background  beam-plasma  is  given  by 

4U  =  /d3?/  dt  E(?,t).}(?,t)  ,  (81) 


where  the  Fourier  transform  of  the  field  is  given  by  equation  (21  )  in  an 
earlier  report,4  namely 


Eki  *  ~ “  Gij<k>jkj  •  (82) 

on  +  i6  J  J 

In  equation  (82)  G—  is  the  linear  photon  Green's  function  of  the  background 
beam-plasma.  Rewriting  equation  (81)  in  terms  of  the  Fourier  transformed 
quantities  one  obtains 


4b l  =  <2*>4  J  Eki.j_ki  .  (83) 

Next  substituting  equation  (82)  in  equation  (83)  prod  ices 

eabs  =  ~i(2w)4  /  —  Gi j (k ) jkj j-ki  •  (84) 

If  the  current  J  is  to  be  real,  then 

3-k  *  3k  •  (85) 


1A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequil ibrium 
Plasma,  Fiz .  Plazmy,  1_  (1975),  673  [Sov.  J.  Plasma  Phys.  1_  (1975),  371], 

4H.  E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  In  a  Nonequi¬ 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories,  HDL-TR- 
1996  (September  1983). 


24 


Substituting  equation  (85)  in  equation  (84)  results  in 

(0)  ,  ,  /Gi-i  (k  )\  * 

4bs  =  -i(2rr)4  /  dk  ( - -jjkjiki  •  <86> 

\to  +  16/ 

Also,  clearly,  since  the  energy  must  be  real,  then 


4  gr  =  4  8i 


(87) 


Then  using  equation  (87),  we  find 


(0)  _  1  /  ( 0)  (0)4 

eabs  -  2  \eabs  +  eabs  J 


But  using  equation  (86) 

(0)* 


dk 


eabs  ~  /  .  Gi j  ) ikj 3ki 

to  -  16 


Or,  if  we  rename  dummy  indices  {i , j }  as  { j , i  } ,  equation  (89)  becomes 


JO)* 

eabs 


/-G- ji(k)\*  * 

=  -i(2tr)4  /  dk  y~-~rhj  jkjjki  . 

Then  substituting  equations  (86)  and  (90)  in  equation  (88),  we  get 

421  =  -  4  (2n)4  /  dki 


r/Gij(k)\ 

/Gji(k)\  *1 

* 

_Voj  +  i<5/ 

\(o  +  i6/ 

Dkijkj 

the  background  is  to 

be  vanishing 

( 91 )  must 

be  vanishing. 

and  therefore 

3ij(k) 

/Gpi(k)\  * 

to  +  i6 

<o 

+ 

3 

(88) 


(89) 


(90) 


(91  ) 


(92) 


in  which  case  equation  (77)  immediately  follows.  It  is  to  be  stressed  that 
equation  (77)  is  at  best  an  approximation  since  the  background  will  not  be 
completely  linearly  nondissipative.  In  summary  then,  equation  (75)  is  in 
essential  agreement  with  equation  (24)  of  Akopyan  and  Tsytovich*  except  for  a 
typographical  error  of  omission  of  a  factor  of  -2ea  in  the  latter;  also 
Akopyan  and  Tsytovich  there  implicitly  assume  principal  values  with  respect  to 
single-wave  particle  resonance  denominators  as  well  as  the  approximate 
nondissipation  relation,  equation  (77).  However  it  is  to  be  stressed  that 
equation  (75)  is  independent  of  these  assumptions. 
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Proceeding  with  the  further  reduction  of  equation  (75),  one  first  defines 


|  ^m£ij 

\ 

1 


(k1 ,-k,k,k1  )  -  B  |  es [Zmlij (kl ) 


+  Emlii^k1 '~k'ki'k) ]  • 


(93) 


Again  the  B  operation  selects  only  the  part  leading  to  the  bremsstrahlung  form 
given  by  equation  (11).  Using  equation  (93)  in  equation  (75)  produces 


_*arad( la)  _  (  ~ea\  ,  ^  dk1  k1  iro(0)i2 

a  \i6tt3  /  1  co2(Wl  +  i6)4  '  k  ' 

x  votpGpm(k1  )vanGjn(k1  )ek£eki  ^ (U1  -  ^1  *va)  <94> 

*  ^mAij (ki '"k'k'ki  )  • 

If  we  substitute  equation  (73)  in  equation  (22),  then  equation  (93)  becomes 


^n£ij (k1 '~k'k,k1  ) 

-  B(-i )  l  e%  J 


d3Ps 


vsm 


(  2 tv)  d  (jui  -  ki  *vs  +  i  6 


x  [(-a)  +  k*vs]VpsJl  +  vs£(-k-VpJ] 


(  ?  \  T 

L  w  +  0)i  -  [k  +  k-]  J  *vs 


-  iTr6(a)  +  to-,  -  (it  +  it-,  )*vs)  j[(u)  -  ft*vs)Vpsi  +  vsi(k*^Ps)] 

*  k  ♦ — -vPs  )1 

L  oil  -  ki  »vs  +  16  s  J 

+  [(a)!  -  ici*vs)Vpsj  +  Vgj  (jti  *^ps)] 

x  rVp .  +  — ^ — kMpl|f«(o)  . 

L  Psi  a)  -  k*vs  +  i6  PsJ)  Ps 


(95) 


If  we  compare  equations  (94)  and  (95)  with  equation  (11),  it  is  evident  that 
only  the  part  of  equation  (95)  proportional  to  6(o)  +  a>i  -  (it  +  k!)*vg)  can 
lead  to  the  bremsstrahlung  form  given  by  equation  (11),  which  is  proportional 
to  6(o)  -  <*va  +  (<  -  )t)«vg).  Note  that^because  of  the  factor^ 6(u)^  -  >c -)  »va)  in 
equation  (94)  one  can  set  =  k!  *v^  in  6(a)  +  u>i  -  (k  +  ki  )»v^)  and  also 
rename  the  integration  variable  )t1  as  -<  to  obtain  S(a>  -  K»va  +  (k  -  l?)»vg). 
Therefore  the  B  operation  in  equation  (95)  selects,  out  of  the  sum  over  s, 
only  the  term  s  =  0,  and  also  only  that  part  proportional  to  6(oi  +  o>1 
-  (ft  +  )*vg).  Therefore  equation  (95)  becomes 


^mUj  =  B(_i)eB  /  - rr  [(-"  +  **''b)vp3£ 

m*i:)  p  (  2tt )  3  U)1  -  k-j  «vg  +  1 6  px' 

+  vgA(-^*^pg)](_ilT)Di3  (k'k1  -  £e)fPB0)  » 

where  the  operator  (k,k-j  ,j$3^  is  defined  by 

^i j (k»ki r  Pg ) 

=  6(oj  +  ag  -  ()?  +  )?i  )*vg)|[(<o  -  ic*vg)7pgi 

*  ''ei(*,5pB)]  7Pej  *  ^TiTT^TTs  k*’‘?ps] 

♦  K«i  -ir»e)vP8j  +  vsj  (*i  *^pg)]  j^pgi  *  „  _  f.'JjTTs 

Simplifying  equations  (96)  and  (97)  produces 

=  _  u  r  d3Pe  vBm  __  rr...  _  i+.y* 


Em  9  i  35  Eire 


Hire g  /  ^ - -  [(»  -  j?-vB)VPM 

°  (2n)3  -  ki»vg  +  i<5  px 

+  v3£k*^pg]DijfR(g°) 


1UU 

Dii  =  6fu  +  u>i  -  (it  +  k!  )*vp){(«  -  k-v6)  -grA—  +  (<*>1  "  M -v£ 


+  (“  -  ^*^^vP3i  ~  .  ^ 4”+"T5  +  vei(k*vPe)vP3j 


+  vBi(k4p  ) 


(*i  ^pj  +  (“i  -  ^i*ve)vp6i 


+  vBUK*vPgJ  — ^  **»>  T  ^  0)  -  k.v3  +  16 

x  (k-VpB)  +  vgj^^pgJVpg.  +  vpjfo^pg)  -  _  ~  —  Mpe)j  ' 

where  the  arguments  of  D^j(k,ki,p3)  are  suppressed.  If^we^use  the  property  of 
the  delta  function  in  equation  (99)  to  replace  u>i  -  )?i  *v3  by  -io  +  k*v3,  the 
first  two  terms  cancel  and  equation  (99)  becomes 


Dij  =  <$(“  +  “1  “  (it  +  itl  )*V0)j^(io  -  k*vg)vpgi 
+  vgi(k*$p  )vp  +  vgi(k*^p  )  - 


7T  ki  *7e 


*  *  rr  >pp 

o)1  -  ki  •  vg  +  16  p 


VVP«  ♦  vSi(^^p3)  -TIi  V  7  -  (Mp*) 


(100) 


-  (a)  -  k*vg)Vp 


pej  0)  -  k.v«  +  is  L  pe 


*  vBj(J',?ps)7Psi  *  vBj(M-5PS)  --rr:;e ;  ^  (J,W  • 

Using  the  properties  of  the  derivative  operator  and  using  the  delta  function, 
one  may  rewrite  equation  (100)  as 


(101  ) 


Djl j  =  <*(“  +  «l  -  (it  +  iti)*vp)|(u  -  k*  vg)  Vp  [ - qr-^ - -)  ki *$pfl 

I  P1\a)1  -  ki  *vg  +  i 6/ 

-  vSJ(*V?P6)vpsi  *  v6i(J-?Ps)Vpej  ♦  _  lBhs  t  i{)  t,  -*pe 

(101  ) 

-  (<o  -  k-vg)  VP6j(^TI^TTl)'k4p3  "  vei^*^PB^VPBj 

+  v8j(ki*^pg)vpgi +  ^(kT^)  ~  :  ^.+^+  i6  (k^P3)|  • 

In  equation  (101)  the  second  and  seventh  terms  cancel  and  the  third  and  sixth 
terms  cancel.  Further  simplifying  equation  (101)  produces 


^P3)7P3j 


Dij  =  $(<•)  +  «i 


-  (it  *  it,).;e)L  -  iE.5e)k  (. — pi* — -Yk.?t 

B1\ui  -  ki  •  vg  +  i 6J  * 


k«vg  +  i6 


ki •va  +  io 


ffi 


using  the  delta  function  to  replace  o)-|  -  k^  *vg  by  -co  +  k*vg  in  the  denomina 
tor,  we  can  rewrite  equation  (107)  as 

Dij  =  6(w  +  ^  -  (k  +  k-|  ) »vg ) | - - - V-* - — T  (& ij - -  k*vg) 

Wgmg(-a)  +  k  »vg  +  1 6  J  l_y  c2  J 


*1  ‘Vgvgi 


(k.vg)^  *vg) 


to  +  k  «v 


When  terms  are  combined,  equation  (108)  becomes 
Dij  =  6(oi  +  a),  -  (it  +  k,  )-vg)  (-  (fiij  -  — 


(^•vg)(k+1  .$.) 


oj  +  k  *v 


u)  +  k  »v 


o)+k»VQ  +  i6 


a)  +  k  »v 


(w  -  k»vg  +  i 6 ) 


Further  simplifying  equation  (109)  produces 


Dj_ j  =  6(u>  +  w,  -  (k  +  ki)«vg)  1 


v  3i  v  3  j 


v3iv3j 


13  c2 


— -  J C 00  ■  k  *vg)2 


v  bl v  PT  >  r  +  +  \ 

v3jk1i  “  vSikj  +  - -  vS*(k  "  kl^ 


(to  -  k*vg) 


+  v0iv8j  lk1  *k 


k  *v gk i  #vg' 


-► 

k 


|_(u)  -  k*vg  +  i6)2  (a)  -  k*vg  -  i 6 ) 2J 

Equation  (110)  may  be  rewritten  as 


+  + 


pe  • 


Dij  =  -6(«  +  0)!  -  (k  +  k,  )*vg)  y^Tg  |(«  -  k  •vg)26ij 
+  (“  -  ^eHvgikj  -  vgjku) 

- [(«  -  k  *vg)2  +  vg*  (k  -  )(u>  -  k*vg) 

cz 

«  -*■  +  +-►  ^  k 

+  k'kic2  -  k-vgk! -vgjl  - -  ^  +  -  -  . 

J  ^0)  -  k*vg  +  lOj*-  [id  -  k«vg  -  i6Jz 

Further  simplifying  equation  (111)  produces 
Di:j  =  -6(oj  +  oi!  -  (k  +  k,  )-vg)  |(u>  -  k»vg)26ij 

+  (“  ~  k*vg)(vgikj  -  v g j k i i  ]  -  V  -^2-gj  [k*k-,c2 


+ 
ki 
--^r-Tr 


+  w(u>  -  k,.vg-  k*vg)]|  - - +  — - - —  - 

I  -  k*vg  +  l  6  [to  -  k»vg  -  l  6  Jz 

Using  the  delta  function  to  replace  (u>  -  )?i  »vg  -  )?*vg)  by  (-u).|  ) 


(112),  we  obtain 

ij 


Dij  =  -«(u  +«!-(£+£,  )*Vg)  y^ 


(“  -  k,vg)2<5ij 


+  (“  -  k*vg)(vgikj  -  vgjku)  -  vgivgj^k.k,  -  !~j 


-► 

k 

7  ?  'F  . 

[oj  -  k»vg  +  i 


k1 


i  6 ) 2  (a)  -  k*vg  -  i 6 ) 2 J 


PB 


(110) 


(111  ) 

W  ' 

(112) 

w  • 

equation 

(113) 


Substituting  equation  (4)  in  equation  (113)  produces 
Dij  -  -  “  «(“  +  <*>1  -  (k  +  k1  )*V0) 

x  [Aij  (k,  ,k  X  +  AiJ)(ki.k)lt1]*VpB  . 

Next  substituting  equation  (114)  in  equation  (98)  results  in 


+  vgJl)t*Vp0]6(aj  +0)-,  -(£+£,  )*vg)  (115) 

x  (ki.k)*,]-«PBf»«601  . 

Integrating  equation  (115)  by  parts,  then  (since  the  background  distribution 
fpjj°>  is  vanishing  at  infinity),  one  obtains 


Zmtti  =  Bire3  / 


d3P8  j_v  (oi  -  £-vg)vgm 

(2ti)3  +  k-|  *vg  -  i 6 

v3mv0I.kn  t(  .  r,+  ?  \  \ 

- T~: * - ^7  +  “1  ~  (k  +  k1  J*veJ 

— +  k-|  *vg  -  16 


(116) 


X  [4f»*(kl,k)lt*  a{|' (»„»)«,  J.tjgi0’  • 


"  L“X]  ■  i»1;j  ri'^nj  'p^pg 

Performing  the  derivatives  in  equation  (116)  using  equation  (105),  and  then 
using  the  delta  function  to  replace  o>i  -  j?i  *vg  by  —to  +  J?»vg  in  the  denomina¬ 
tors,  we  obtain 


vBnvB£ 


vBnvB£ 


vBnvB£ 


vBmvB£knk1] 


(J?  +  It-,  )«$g)[Aj np;"  (k,  ,k)£  +  A{j  (k,  ,k)lt!  ]»7p  f 


V -.-vwvv v\.-- 

\  «/“  »'  ,* 


Combining  terms  in  equation  (117)  produces 


k  *vBk1  *vg 


to  -  k*vg  -  i6  (to  -  k*vg  -  ifi)2  (to  -  k*vg  -  i 5 ) 2 
x  6(oj  +«!-(!?  +  (ki,k)f?  +  (ki  ,k  )#i  ]  *^pgfpg0) 

Further  simplifying  equation  (118),  we  obtain 

-  ,  r  d3Pg  1  x  kmv3Jl  “  kUv0m 

^mfij  -  J  (2ir)3  Tptrg  \-S*i  ~  ~  _  ^  _  i6 

+  T? - VfmIg£  — -rj  [<*>(<*>  -  (k  +  k,  )*ve)  +  k*klC2] 

c^(to  -  k*vg  -  i6Jz 

x  6(d)  +  u>i  -  +  *vg)  [a{  j  *  (ki  ,k  ))t  (119 

+  Aij,(k1,k)jr1]^ppfJ<e0)  . 

If  we  use  the  delta  function  in  equation  (119)  to  replace  (to  -  ()?  +  i?-|  )*vg) 
by  (-u^ ),  then  after  rearranging  terms,  we  obtain 


^m£ii  ~  ~Bire«  / 


0  (2ir) 


vB*v3m 

7 - T~' * - 

[to  -  k»vg  - 


PB  1  r  vB*km  -  v&nk  1  £ 

r?  ^  L  * » -  iUs  - 1« 

77TF  (**'  -  9)}(“ +  -  -  <*♦*>>•*•> 


X  [Aif)#(ki,k)jr  +  illj,(ki,k)<1]-Lf 


PB  P0 


Next  substituting  equation  (4)  in  equation  (120)  produces 

imiij  -  -01ieg  /  Alta)(k1'k)6(“  +  “1  -  (k  +  *1  Wp) 

X  [Aif*(k,,kX  +  Aiftk^k)!?!  ]^Pgfpg0>  . 


1 


■••.V.  ^vSVv.vv  -v;. 


.  S  .-  ■.»' 

*  •  •  *  *  -  ’  • 


i<S 

<  A  I 


Equation  (121)  is  to  be  compared  with  equation  (26)  of  Akopyan  and  Tsytovich. 1 
Since  the  B  operation  projects  out  s  =  8  as  discussed  above,  equation  (93)  may 
be  written 


^mJlij  (^i  »-k,k,ki  )  -  Beg[zm^j  (k-)  )  +  Zm£ji  (ki  ,-k,ki  ,k  )  ] 

Comparing  equations  (121)  and  (122),  we  find 
B[s»f!j(*1»-k.k.k1)  +  4g]i(ki,-k,ki,k)] 


(122) 


=  -Bireg  / 


6  J  (2ir)  3 


S(oj  +oji  -  (k  +  ki  )*vg) 


(123) 


*  *  a!3si (*,.*)>?,  Hpgf^0’  . 


By  equation  (4)  it  follows  that2 


6(io  +  0)!  -  (j?  +  £•,  j-vgjA^5  (k,  ,k )  =  6(o>  +  tu1  -  (jt  +  ^  )»vg)A^f>*  (k.k,  )  . 

(124) 

Substituting  equations  (76)  and  (124)  in  equation  (123),  dividing  by  i,  and 
taking  the  complex  conjugate,  one  obtains 


(125) 


®  [^mfi  j  (ki  ,k,  -k,  -ki  )  +  T^  (k,  ,k,-ki  ,-k)  ] 

=  Biireg  /  — ~  6(u  +  oil  -  (k  +  ^  )*vg)A^)  (k,ki  ) [a^  (kt  ,k )k 

(2TT)3 

*  Aif*(k,.*)e,]-5P64(S0)  . 


If  one  ignores  single-wave  particle  resonances,  the  i6  in  equation  (4)  may  be 


ignored  and  A>P'(k1fk)  becomes  effectively  real  and 

plAt?,*(ki.k)  =  P1A^)(k1,k)  , 


(126) 


where  denotes  the  operation  of  taking  principal  values  with  respect  to 
Cerenkov  resonance  denominators  of  the  form  (oo  -  Jt»Vg  +  i 6 ) .  The  principal 
part  of  A^^(kj,k)  is  real.  The  nonprincipal  part  is  proportional  to 
6(<u  -  £*Vg]  or  6'(o)  -  £*Vg)  according  to  equations  (4)  and  (73).  Therefore 
it  will  not  yield  the  correct  bremsstrahlung  form  and  may  be  dropped. 
Operating  on  both  sides  of  equation  (125)  with  P1  and  using  equation  (126),  we 
obtain 

V.  Akopyan  and  V .  N .  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plazmy,  1_  (1975),  673  [Sov.  J.  Plasma  Phys.  1_  (1975),  371], 

2H.  E,  Brandt,  Nonlinear  Force  on  an  Unpolarized  Relativistic  Test  Parti¬ 
cle  to  Second-Order  in  the  Total  Force  in  a  Nonequilibrium  Beam-Plasma  System, 
Harry  Diamond  Laboratories,  HDL-TR-1995  (August  1983). 


P1  B[Tmfci  j  0*1  )  +  Tmfcji  0*1  •  **'  1  •  )  ] 

d  3  *► 

=  PiBiireg  /  — —  S(w  +  to-j  -  (it  +  ki)*vp)  (127) 

(  2  it  ) 

-  aJI’M.  >!?’(*, .*)(?  *  . 

Equation  (127)  agrees  with  equation  (26)  of  Akopyan  and  Tsytovich1  except  for 
k  appearing  interchanged  with  k1  in  and  A^j^.  This  is  apparently  an 

error  in  the  cited  work.1  The  operation  P^B  is  implicit  there  and  the  i  is 
absent  since  the  imaginary  part  has  been  taken  (the  coefficient  of  i)  and  the 
right-hand  side  of  equation  (127)  here  is  purely  imaginary. 

Next  substituting  equation  (121)  in  equation  (94)  results  in 


+crad(la) 

<Fo 


e«e|  ,  d3pg  dk  dk-| 

16tt2  ( 2ir ) 3  uj2  ( cu-|  +  ifi)1* 


x  jc-)  Ej^101|26(u3  +  -  (it  +  ici  )*vg)6(u1  -  it-,  *va) 


x  votpGpm(k1  )A£m  (ki  »k)eki[ekiAij  (k  1  »k)Gjn(k1  )vank 

+  ekiAij)(k1*k)Gjn(kl)van^l]*^pBfpg°)  • 


(128) 


Using  equations  (92)  and  (124)  in  equation  (128)  produces 
<forad(1a)> 

.  B  fM  r  d3^  —  dkl  i  Ieo(0)|2 


1  6ir2  (2n)3  u)2(o)1  +  i6)2tti2 

x  6(u  +  <*>!-  (it  +  )-v6)6(uj1  -  lti«va)[e^£AI{1f)(k,ki)Gmp(k1)vap]*  (129) 


x  [^iA5i,(k,k1)Gjn(k1)vcmlt  +  e^iA^)*(k,k1)Gjn(kl)van^l]^p3f?g0)  • 
Integrating  over  ,  letting  iti  =  -it,  and  denoting 

<  =  (*•$„,*)  , 


(130) 


then  we  find  that  equation  (129)  becomes 


V,  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plazmy,  ±(1975),  673  [Sov.  J.  Plasma  Phys.  ±  (1975),  371). 


O.V,V*,;v 


••  •••  .*•  /■  -V  . 


*.  •*-**»  '* .  *' .  •*«  v »*.  •*. 

•  *#«  •.»  *  *mn 


-B 


♦orad(la) 
<Fa  > 


e2e2 

eae6 


16* 


p  d3pg  du  d3k  d3ic  k  E^°^|2 
2  (2*)3  u)2(K»va  -  i6)2(tc»va)2 


x  6(a)  -  +  (<  -  it)*ve)[eg£A1ilf)(k,-(c)Gmp(-K)vap]* 

x  [eki^ji  (k, -<)Gjn( -<)vank 


(131  ) 


> 

•V 


fR(0) 

P3rPg 


If  Cerenkov  resonance  is  ignored,  equation  (126)  applies.  This  is  justified 
since  only  the  principal  part  of  (k,-ic)  leads  to  the  bremsstrahlung  form, 

and  the  principal  part  is  real.  JAlso,  ignoring  the  i6  in  (ic»va  -  i6)“3, 
dropping  the  unneeded  B  operation,  and  making  the  change  of  indices  from 
{i,j,n,p}  to  {j,i,p,n},  we  find  that  equation  (131)  becomes 

<+arad(1a)>  =  e°eg  ,  d3Pg  d“  d3*  d3*  ^ | Ek ( ° 5 ) 2 
16*  2  ( 2ir )  3  oj2(  <•  Vq)1* 

x  6(a)  -  k#Vq  +  (k  -  it)»vg)  (132) 

x  |ekjAij  )(k,-K)Gip(-K)vapj2(ic  -  ic)*Vpgf^0)  . 

Equation  (132)  is  in  fact  equation  (3),  as  was  to  be  shown.  It  is  to  be 
stressed  that  to  obtain  this  expression  it  was  necessary  to  assume  that  the 
background  is  approximately  nonabsorptive  and  spatially  isotropic.  However, 
these  assumptions  are  not  involved  in  equation  (128).  As  already  stated  in 
the  introduction,  equation  (27)  of  Akopyan  and  Tsytovich3  is  apparently  in 
error;  namely,  a  factor  of  2  is  omitted  there  and  the  arguments 
of  A>j'(-K,k)  are  interchanged.  The  additional  factors  of  (4*)2(2*)"6  there 
are  apparently  due  to  the  use  of  Gaussian  units  and  different  Fourier 
transform  conventions,  as  already  discussed  in  equation  (75).  Also  the 
Green’s  function  used  there  is  symmetric. 
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3.  CONCLUSION 

A  number  of  important  physical  assumptions  and  mathematical  techniques 
have  been  identified  which  are  involved  in  Tsytovich’s  theory  of  nonlinear 
bremsstrahlung  and  radiative  instability  in  relativistic  nonequilibrium  beam- 
plasma  systems.  In  particular,  the  stochastic  average  of  the  lowest-order 
Born  approximation  for  the  nonlinear  force  on  the  bare  charge  of  a  relativ¬ 
istic  test  particle,  equation  (1),  has  been  reduced  to  the  form  given  by 
equations  (2),  (3),  (47),  (25),  and  (26).  In  the  result  of  Akopyan  and 

Tsytovich,1  their  equation  (27)  which  corresponds  to  equation  (3)  here,  a 
factor  of  2  is  evidently  omitted  and  the  arguments  of  are  incorrectly 

interchanged.  Equation  (3)  has  been  shown  to  hold  providecT  the  background  is 
approximately  nonabsorptive  and  spatially  isotropic.  However,  the  form  given 
by  equation  (128)  holds  without  these  assumptions.  Also,  an  important 
expression,  equation  (123),  was  obtained  involving  the  third-order  nonlinear 
conductivity  tensor.  Using  this  expression  it  was  shown  that  in  equation  (26) 
of  Akopyan  and  Tsytovich1  the  arguments  of  are  again  incorrectly 

interchanged. 

The  present  work  will  facilitate  better  understanding  of  the  methods 
necessary  to  reduce  the  other  components  of  the  total  nonlinear  force  due  to 
bremsstrahlung  emission. 1  The  latter  is  needed  to  determine  the  nonlinear 
bremsstrahlung  probability  and  the  conditions  for  the  occurrence  of  a  brems¬ 
strahlung  radiative  instability  in  nonequilibrium  relativistic  beam-plasma 
systems. 


lA.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plazmy,  1  (1975),  673  (Sov.  J.  Plasma  Phys.  1  (1975),  371). 
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